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In this paper the interaction of high energy neutrinos with weak gravitational fields is briefly 
explored. The form of the graviton-neutrino vertex is motivated from Lorentz and gauge invariance 
and the non-relativistic interpretations of the neutrino gravitational form factors are obtained. We 
comment on the renormalization conditions, the preservation of the weak equivalence principle and 
the definition of the neutrino mass radius. We associate the neutrino gravitational form factors with 
specific angular momentum states. Based on Feynman diagrams, spin-statistics, CP invariance and 
symmetries of the angular momentum states in the neutrino-graviton vertex, we deduce differences 
between the Majorana and Dirac cases. It is then proved that in spite of the theoretical differences 
between the two cases, as far as experiments are considered, they would be virtually indistinguishable 
for any space-time geometry satisfying the weak field condition. We then calculate the transition 
gravitational form factors for the neutrino by evaluating the relevant Feynman diagrams at 1-loop 
and estimate a neutrino transition mass radius. The form factor is seen to depend on the momentum 
transfer very weakly. It is also seen that the neutrino transition mass radius is smaller than the 
typical neutrino charge radius by a couple of orders of magnitude. 

PACS numbers: 13.15.-|-g, 04.25.Nx, 14.60.Lm 



I. INTRODUCTION 

Neutrinos are very weakly interacting particles, experi- 
encing only the weak nuclear force and the gravitational 
force. They arc produced copiously in many high energy 
astrophysical processes and due to their weakly inter- 
acting nature travel almost unhindered to earth. These 
characteristics make them a potentially powerful probe 
of astrophysical phenomena which are otherwise inacces- 
sible. It is speculated that these ultra high energy neu- 
trinos may have energies as high as ~ lO^^GeV. (See for 
example [l| and references therein.) Many ongoing and 
planned experiments aim to detect these high energy as- 
trophysical neutrinos 0- Results from these neutrino 
telescopes have the potential to provide deep insights in 
cosmology and particle physics. 

Our aim is to study how high energy neutrinos, with 
3> rrii,, interact with gravitational fields that are 
weak. Many of the gravitational systems that we arc 
generally interested in have weak gravity. To our knowl- 
edge the first study on neutrino gravity interactions was 
that of Brill and Wheeler Q in the context of introducing 
spinors in general relativity. There have been subsequent 
studies on neutrino-gravitational effects, mainly pertain- 
ing to neutrino oscillations. 



'Electronic address: I aamenon@umich.edu I 
t Electronic address: |madhav@uchicago.edu| 



One approach pioneered in [J| was to assume that 
the weak equivalence principle is slightly violated in na- 
ture, meaning that the gravitational coupling constant 
is non-universal and depends on the flavor of the par- 
ticle species, and then probe its effect on oscillations. 
Another approach has been to dcflnc a covariant gravita- 
tional phase Q (in the context of the Dirac equation in 
curved space-time) that would be the relevant substitute 
for the vacuum phase e*^^ in curved space-times and ex- 
plore its effects on neutrino oscillations. There have also 
been detailed studies of gravitational effects on neutrino 
oscillations in a medium [6[ and on spin oscillations in an 
external gravitational field ■ 

Towards the goal of understanding neutrino gravity 
interactions we first study the graviton-neutrino vertex 
using general symmetry principles and then explore the 
non-relativistic limits of the neutrino-graviton form fac- 
tors. This line of reasoning was motivated by the use of 
invariance principles, long back, in studying the photon- 
neutrino coupling 0]. We use arguments from Feynman 
diagrams, Fermi-Dirac statistics, CP invariance, cross- 
ing symmetry and symmetries of the angular momentum 
states to study the vertex. This would, among other 
things, give us some insight into how the Majorana neu- 
trino and Dirac neutrino cases could be different as far 
as gravitational interactions are concerned. Remarks are 
also made on the weak equivalence principle and its valid- 
ity in the context of renormalization. It must be pointed 
out that some of the differences between Majorana and 
Dirac form factors were anticipated earlier in [oj using 
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methods different from the one we follow. 

In spite of theoretical differences, we prove to 
0{m,y/Ei,) that the Majorana and Dirac cases cannot 
be distinguished by gravitational interactions, as far as 
experiments are involved. This result is found to hold 
true for any space-time geometry that satisfies the weak 
field criterion and is a simple extension of the practical 
Majorana-Dirac Confusion theorems Q first discussed in 
the context of neutrino scattering and neutrino-photon 
interactions. 

We also perform an approximate calculation of the off- 
diagonal (in mass basis) neutrino-graviton form factors 
at 1-loop to understand their dependence. Using the 
loop calculation results the Vi Vj neutrino transition 
mass radius is estimated when i ^ j. A comparison is 
then made of the neutrino transition mass radius and the 
neutrino transition charge radius. Finally, the difference 
in order of magnitudes of the two radii is motivated based 
on some physical arguments. 



II. THE NEUTRINO-GRAVITON VERTEX 



For weak gravitational fields we may expand the metric 



as 



(1) 



where k = V327rG, rj^y = (1, —1, —1, —1) and h^^" is in- 
terpreted as the spin-2 graviton. G is the Newton's grav- 
itational constant. The usual criterion for considering a 
gravitational field to be weak is 5g ~ <C 1 

For the case of spherical symmetry, the exterior is the 
Schwarzschild spacetime. In this particular case 

2GM 



c^R 



where M is the mass of the body, R is the radial distance 
and c is the speed of light. Let us consider some typical 
examples now. For the exterior of the sun the above 
expression gives 
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Active galactic nuclei typically have M 
i? ^ 10 KPc which means 



IQI^Mq and 
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A typical neutron star mass and radius, range between 
1.3 - 2.1 A/q and 10 - 20 Km [13]. This typically leads 
to 



0.3 



Thus we may be optimistic that the analysis we perform 
is valid to a good extent in many astrophysical cases of 
interest. 

The graviton-neutrino vertex current J7^i/ in Fig. ([T|) 
transforms as a symmetric rank-2 tensor. The curly line 




FIG. 1: Neutrino Graviton Vertex. 

denotes the off-shell graviton. The spin-2 gauge current 
is just the energy-momentum tensor of the neutrino in 
the gravitational background. By Lorentz invariance the 
possible vertex factors are 

u{p )iij^^u{p)i = u{p )i A{q^)gi_,^ + B{q^){ri_,r„) 
+C{q^){q^q,) + b{q^){^{^r,}) + E{q^){^{^q,}) 
+F{q^){q{^r,})\u{p), (2) 

where p is the incoming neutrino momenta, p is the out- 
going neutrino momenta, q = p — p, r = p + p and 
{} denotes complete symmetrization of the indices. The 
subscript i in the spinors label the mass eigenstate of the 
incoming and outgoing neutrinos. 

Using the Slavnov- Taylor- Ward identities 

qf^ u[p )i ij^^ u{p)i = 

to impose gauge invariance in Eq. ^ we find that the 
vertex in general must have the form 

u{p )t'iJ^,„ [i, i] u{p)i = u{p )j Fi{q^){q^g^^ - q^,q^) 
+F2{q^){r^r,) + F^{q^){^{^r,})\u{p), . (3) 

To arrive at the above expression we used the Dirac 
equation along with the result r -g = 0. The notation [i, i] 
denotes the mass diagonal case. Also, we will label the 
parity conserving and parity violating currents by J^^^, 

and J^^, respectively. The form factors we have defined 
above may be shown to be linear combinations of the 
ones defined, for the mass-diagonal case, in |ll| in the 
context of studying the gravitational decay of a spin-1/2 
particle. 

It will be seen that when the relevant graviton Feyn- 
man rules are derived to lowest order in k, only the 
7{q''/3} operator appears at tree level. Also note that 
the operator r^ri, is a chirality flipping operator. Due 
to the tiny neutrino masses and the left handed nature 
of the electroweak currents, this would imply that this 
operator is suppressed by mu/ E^. 
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Until this point we were considering only parity con- 
serving operators. If we include operators with 75 and 
proceed as above, we find that there are again three pos- 
sible vertex operators 



+G2{q^){r^r,)j' + G3{q'){qH^^^r,}) 
-2m^.{q{^r^})}-f^ u{p), . 

For the off-diagonal (in mass eigenstates) case ui 



(4) 



the only non-zero component of J7^^ is jToo and the above 
expression reduces to 



K h°°u{p^ 0)1 iJoo u{p^ 0)i ~ K h^°u{p-* 0), 
4(F2(0) + G2iOh^)ml^ + AF3iO)mJ u{p ^ 0), 



Based on the above expressions let us re-define the form 
factors with appropriate normalization factors so as to 
obtain a simpler result. 



with mass eigenstate i ^ j similar arguments immedi- 
ately give 

Hp ["iJ] Hp)i - Hp )j Ei{q^){q^g^,^ - q^q^) 

+E2{q^){ry,r^) + E^{q^){q^{-i{^r^}) - Ay mi,(gf{^r^})} 
u{p)^ + 0{^^Jml) , 



en / 2n ^2(9^) „ , 2n G2{q^) . 2x 

^2(9 ) -r—, G2{q ) ^ -r—, F^{q ) 



Fsiq') 



4m „ 



4m „ 



(5) 



Hp )jiJ,,^ [hj]Hp)'i - Hp )j ^1(9^)7 



In the above expressions, it is assumed that m,^. are neu- 
trino masses measured from some kinematics. Thus they 
may be regarded as neutrino masses in the 'inertial' sense. 
With this new normalization we then have 

K/i°°i2(0),iJoow(0), « m,,$, [F2(0) + F3(0)] 0^ (7) 

where (j) is the usual "large component" Pauli spinor (of 
the neutrino) and <i>j, is the gravitational potential in the 
neutrino rest frame. To arrive at the above expression we 



{q g^iiy - Q^iQi^) + D2{q ){r^r„)j + Daiq ){q {ji^r,,}) used the well known fact that in the weak gravity limit 



(6) 



sitions. Also, we define Aj 
rriv, +miy, and A^m?, = m?,. — m 



- T.,jm^{q[f,r^})}j u{p)i + 0{A,-jm^) , 

where the notation denotes mass off-diagonal tran- 

I • writing the above 
expressions for the off-diagonal case we have neglected 
terms of ©(Ay-mJ). Since we are not interested in any 
neutrino oscillation phenomena this is justified. For the 
exact vertex in the case oi i ^ j the reader is referred to 

Let us now try to deduce physical interpretations for 
the neutrino form factors in Eqs. Q and ^ in the non- 
relativistic limit. In the electromagnetic case it is well 
known that in the non-relativistic limit the form fac- 
tors are interpreted as charge and electromagnetic mo- 
ments. For a neutrino in that case the neutrino charge 
vanishes. It is thus interesting to see what the interpre- 
tations would be for the neutrino-graviton case. It must 
be mentioned that the energy-momentum structure form 
factors and their non-relativistic interpretations were first 
studied in the context of understanding the 'mechani- 
cal structure' of particles and the self-energy of the elec- 
tron [3. 

Consider the neutrino form factors in the limit q ^ 0. 
Then we have the graviton-ncutrino interaction Hamilto- 
nian density 

limT^g = Kh'^" Hp)i'>-^tiiyHp)i 

~ Kh'^-u{pU^[F2{Q) + G2{0h')p^.p. + 2^^3(0)7{mP.}] 
u{p), . 

Further simplification is achieved if we transform into the 
rest frame of the incoming neutrino. In the rest frame 
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,00 



1 + 2$, 



in natural units. 

From Eq. ([7]) we deduce that F2{q^) and ^3(5^) are 
associated with the energy density or mass of the incom- 
ing neutrino. This observation leads to the requirement 
that in the limit g — > we must have 



J^2(0)+ F3(0)= 1 



(8) 



But we already noted that r^r^ is a neutrino chirality 
hipping operator. In the q ^ limit ^2(9^) niust vanish, 
since in the limit of zero momentum transfer we do not 
want to induce any neutrino chirality flips. This may 
be considered as a consequence of angular momentum 
conservation. Hence, 



F2(0)^0 



(9) 



The above conditions imply that any correction, to the 
neutrino mass, calculated at one-loop be renormalized 
such that in the limit g it vanishes: 

AF, (<?2),,„„^ = AF,(q2)^^^^^_ Af^,(0)i,„„, 

AG, (g2),,„„^ = AG,(g2),,^„^_ AG,(0)i,„„p . (10) 

This renormalization condition is exactly analogous to 
the electromagnetic case where in the limit g ^ we 
require the one-loop vertex corrections to not renormalize 
the electric charge or moments. 

Note also that for the off-diagonal case when i ^ j 
such a renormalization must again be imposed as q ^ 
to prevent the generation of off-diagonal elements in 
the non-relativistic limit. This condition will be used 
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explicitly later in section IV when we calculate the mass 
ofF-diagonal form factors. 

A version of the weak equivalence principle states that 
the gravitational mass is equal to the inertial mass for all 
particles We see that, imposing the conditions in 

Eq. ([TU)) at every step ensures that the mass eigenstates 
satisfy the weak equivalence principle. This would also 
mean, among other things, that a weak gravitational field 
does not affect neutrino oscillations. The relevant oscil- 
lation in pure gravity (weak) would still be the vacuum 
oscillation. 

Analogous to the electromagnetic case, one can now 
consider the F3{q^) form factor roughly as the Fourier 
transform of a hypothetical neutrino mass density p'^{r), 

This leads naturally to a definition of the neutrino mass 
radius 



(11) 



Later, we will estimate the i ^ j neutrino transition mass 
radius by explicitly evaluating the Fcynman diagrams at 
1-loop. 

Let us now generalize the previous results for the neu- 
trino form factors. By the virial theorem, for any energy- 
momentum tensor we have the general result (see for ex- 
ample [io| ) 



d^'x {J^)t = / d'x {J^), = (E), 



where {E)t is the time- averaged energy and ()t denotes 
time averaging. This implies that 

mt = {PE)t 

is the average energy density. So if we are considering 
neutrino interactions in an environment of gravitons, in 
the non-relativistic limit, we may use the above expres- 
sion in Eqs. ([3|) and ([4]) along with the normalizations 
and interpret 



AF^{q') 



F2{q^) 



4m „ 



To summarize, the complete mass-diagonal graviton- 
neutrino vertex is of the form 

u{p')iijf„yu{p).i = u{p')i Fi{q^){q^g^^ - q^q^,) 
F^iq') ,_^F,{q^)f , 

+Gi{q'){q'g^. - q^^q.h' + ^^(^M^^)7' 
+G3(g'){9'(7{MM) - 2m,,((Z{^r,})}75lu(p), , (13) 



with the interpretations given above. 

Further relations between the form factors may now 
be obtained by including information about the nature 
of the standard model gauge interactions. To be specific, 
if we incorporate the observational fact that electroweak 
currents are left-handed and neutrino masses are very 
tiny then by considering the matrix elements between 
various chiral states we can further constrain the form 
factors. We will see later on that at tree level the only 
operator that appears is j^affsy- The other factors are 
generated at higher orders, so we have the obvious con- 
straints on them from electroweak theory that matrix ele- 
ments of right-right chirality or right-left chirality vanish 
owing to the tiny neutrino masses. For example 







(12) 



would imply (choosing an appropriate non-chiral Dirac 
basis) that 

AFsiq'-) - AG3(g2) 

to Oim^/E^). 

A similar line of reasoning will be used when we address 
what the experimental consequences of the results we ob- 
tain are, for distinguishing Majorana and Dirac neutrinos 
through gravitational interactions. 

Later in our analysis we will also explicitly take the 
non-relativistic limit of the operators in Eq. (|13p and 
identify the angular momentum states associated with 
each one of them. This along with other arguments from 
spin-statistics, CP, crossing symmetry and Feynman di- 
agrams will yield information about how the Majorana 
and Dirac cases differ as far as neutrino gravitational in- 
teractions are concerned. 



as an average energy density for the incoming neutrino. 
The first term is familiar from the (/ — > case and the 
second term seems to give a next order correction in q^. 

All the above expressions give an interpretation for the 
neutrino form factors Fi(g^), F2{q^) and F'g,{q^). The 
form factors Gi[q^), G2{q^) and G3{q'^) may be thought 
of as the 'axial' counterparts of the F{q^) form factors. 



III. ANGULAR MOMENTUM STATES AND 
SYMMETRIES OF THE VERTEX 



Let us now consider the hypothetical case of identical 
Majorana fermions in the final state, for the s-channcl 
process with matrix elements [vij) )iv{p)i\S^v\^) , which 
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is related to the t-channel process (through crossing sym- 
metry) with matrix elements {v{p )i\J^fi,v\v{p)i). Explic- 
itly the process under consideration is 

h*{q) v,{p) + v^{p ) . 

In the above, h* is the off-shell graviton associated with 
the weak gravitational field. For the s-channcl process 
we have the new definitions q = p +p and r = p —p. As 
an aside, note that unlike processes such as ^ v + -f 
where the photon may be on-shell, in a t-channcl process 
Vi + h Vj the spin-2 graviton cannot in general be 
on-shell, due to angular momentum conservation. The 
possibility of an on-shcU graviton considered by K.L. Ng 
in seems to be in disagreement with this reasoning. 

For exploring the differences between graviton- 
Majorana neutrino and graviton-Dirac neutrino form fac- 
tors, we will adopt a line of reasoning similar to the one 
used earlier in the context of neutrino photon couplings. 
In the case of photon couplings, Fermi-Dirac statistics 
and some general results were used to show that a Ma- 
jorana neutrino can have only a single electromagnetic 
form factor compared to four in the Dirac case Q . 

Henceforth, a super-script D will label the Dirac neu- 
trino case and M will label the Majorana neutrino case. 

Assume that the final state fermions in our case are 
identical Majorana neutrinos. By Fermi-Dirac statistics 
we must then have 

{u{p )iv{p)i\J^,u\Q) = Ui(p {p ,s), {p, s) Vi{p) 



Ui{p)iJa,y iP,s)i (P :S ) Vi{p ) 



(14) 



Wc note that under this exchange q q while r — s- —r 
in J^fiv In addition to the above equation wc also have 

uip ,s')ij^'''vip,s) = [i-/2vip)f-/oiJ'"'[-il2U*{p)] 

.1 T 



v{p ) = -u{p) 



v{p ) (15) 



where C ~ 7"7^. For Majorana neutrinos, from Eqs. (fT4] 
and (|15|) . the s-channel matrix elements must satisfy 



)tv{p)t\J^,^\0) = Ut{p ^{p ,s ), (p,s)J v^{p) 

= -u{p)i-iJ^l {P, s), {p ,s) v{p )j 

= -u{p), [c-^^JHI [(p, s'), (p, s)] C] ^ v{p\ . 

Thus wc may conclude from above that for a Majorana 
neutrino mass cigenstate the only form factors that sur- 
vive are those for which 



{P ,s ),(p, s) 



C 



J, 



M 



{Pi s),ip ,s) 



It is to be emphasized that the above conclusion is de- 
rived solely from Fermi-Dirac statistics and the fact that 
for Majorana neutrinos the final state in the equivalent 



s-channel process consists of identical fermions labeled 
by i. Using the above criterion on the terms in Eq. ([13]) 
it is found that for the Majorana case the form factor 



Criq') = 



(17) 



This was also noted previously, following a different 
method, in Q. Thus we may claim that, in the coupling 
of a neutrino to a graviton, the Dirac neutrino has six 
non vanishing form factors while the Majorana neutrino 
has only five form factors to be consistent with Fermi- 
Dirac statistics. 

We also note that for the cases when there are tran- 
sitions ly^^ —> v^^ with i ^ j, by crossing symmetry, we 
have a correspondence between the processes 



A/ h* 



,M 



^ h* 



M h' 



,M 



implying that the form factors of these matrix elements 
are the same. 

A more transparent way to understand why the form 
factor G^^ {q^) vanishes in the Majorana case is by look- 
ing at the angular momentum states (of the final neu- 
trinos) produced by the virtual massless graviton. This 
method also leads to other insights regarding neutrino 
gravitational interactions. 

A virtual massless graviton has a J = as well as a J = 
2 component ■ This is very different from the photon 
case where there is only a J = 1 component, for both on- 
shell and off-shell photons. The possible neutrino final 
states (Majorana or Dirac) for the graviton case are : 



J = 



S = Q : ^So , 



(18) 



S = l 



where as usual J = L-\- S and we have used the standard 
spectroscopic notation ^^'^^Lj for the angular momen- 
tum states. If we take the on-shell limit for the graviton 
some of the states cancel to yield the expected J3 = ± 2 
helicity [l^. Finally, note that the six angular momen- 
tum states of Eq. p8)) are related to the six operators that 
appear in the Dirac neutrino-graviton vertex. This asso- 
ciation will be explicitly derived in the non-relativistic 
limit. 

Let us now consider the Majorana case where the two 
final state particles are identical fermions. We require 
this final state to be anti-symmetric under the exchange 
[pl,sl]i <-> [p2,s2]i. For an identical (Majorana) two- 
particle final state we have under the permutation oper- 
ator Pi2([pl,sl], ^ b2,s2],) : 



^'AC'^'Lj) = (-1) 



L + S+l 



(19) 



Using the above expression for the Majorana final states 
we find that five of the final states are antisymmetric 
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except Z?2 which has 



(20) 



This again imphcs that in the Majorana neutrino case one 
can have only five graviton form factors to be consistent 
with Fermi-Dirac statistics. When we derive the angular 
momentum correspondence of the graviton-vertex oper- 
ators explicitly in the non-relativistic limit, we will show 
concretely that Eq. ([^0]) indeed implies Eq. pT]) . 

Now we will try to associate each of the form factors 
(in the non-relativistic limit) with states of fixed angular 
momentum. The off-shell graviton polarization tensor 
e/ijy [J, J3] with J = and J = 2 may be written in terms 
of the j = 1 , j' = 1 polarization vectors r?^ [j, J3] , 77^ [/ , jg] 
of the spin-1 field (photon) as 

33+j'3 = J3 

(21) 

Here C{J j j 'jJsjaj^) are the usual Clebsch-Gordan co- 
efficients. We will work in the rest frame of the off-shell 
graviton h* so as to simplify the analysis. The photon 
polarization vectors of good angular momentum in this 
special frame may be chosen to be 

1 



r,[l,+l] = -^(0,-l,-*,0) 



?7[1,0] 



(0,0,0,1) 



(22) 



V2 



(0,1, -i,0) 



In this case the five graviton polarization tensors with 
J = 2 may be written explicitly using the above formula 
as 



e[2,±2] 



e[2,±l] 



^[2,0] = :k 



/o o^ 
1 ±i 

±i -1 

Vo oy 
/o 

Tl 

Vo Tl -i 

/O ON 
0-100 

0-10 

Vo 2y 



(23) 



and the J = polarization tensor is 



e[0,0] 



1 



/o 

0-10 
00-10 

Vo -ly 



(24) 



In the rest frame of the off-shell graviton we have the 
neutrino momenta p2 = —pi = p and without loss of gen- 
erality we normalize it to have unit magnitude. Following 



the notation in M we label 



S and (j)2S(j)\ = T for 



the singlet and triplet states. Here and 02 are again 
the large component Pauli spinors of the two neutrinos. 
With these definitions, for the states of interest, we find 
that the angular momentum eigenfunctions take the form 



^^0 



S , 
p-f , 

{3p®p-l®l)S , 
p®f, (p-f) , 
p (g) {f xp) , 



: 5{p-T)p(g)p+{p-T)l®l-p®T. (25) 

Here (g) denotes a direct product of vectors and sym- 
metrization is implicit. Working in a non-chiral Dirac 
basis one can take the non-relativistic limit of the oper- 
ators in Eq. (|13p in the standard way and express them 
in terms of the angular momentum eigenfunctions above. 
Using Eqs. ([2T|) - ([25|l along with the properties of the po- 
larization tensors gives 

e^'iO, J3] ■ uJq^g^. - q^q.)vi ^ -^p ■ f 



NR 2ey 



Hp ■ T)piPj 



10m2 



+ (p- r)(5,j - {P{^TJ}) 

e^''[0, J3] • U2{^i_ir^ + ^^r^vi 

NR 4e'^ r 



4(p.T),5„ + (p{,r,}) 

NR p-T 



e>"'[2,J3]-U2ir^r,)vi 
HP ■ - {p{ifj}) 
e^"'[Q,.h]-U2{r^r,)vi 



NR ^p-T 



5{p-T)p,pj 
{PitTj}) - {p-T)5,j 



(26) 



£^''[0, J3] • U2 {q^Qi^y - qt^qv) l^Vi 



NR. 



e^''[2, J3] • U2 il^r, + l.r^) ^'v, ^ — p{.,{p x f ),} 

Tfli, 



e^'^[2,J3] •U2(r^r,)75t;i 



e^"[0,J3]-U2(rpr,)7V 



NR _4 / j/_ 
' 3 ^ 4to2 

e'-' {ipiPj - dij) S 

NR 



4 1 



S 



(27) 



where, again {} denotes complete symmetrization with 
respect to the relevant indices. 
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Comparing Eqs. ([25|) . (p6)) and ((27|) we deduce the an- 
gular momentum associations, in the non-relativistic 
limit, for the graviton vertex operators 



9m 



(28) 



'F2 



'P, 



'Po 



'Po 



'D2 ; ^So 



This correspondence of the vertex operators with the an- 
gular momentum states in the non-relativistic limit again 
provides an interpretation for the associated form factors 
and is complementary to the earlier results we derived in 
Eqs. and nil). 

Furthermore we see from Eq. pS)) that the state is 
indeed associated with the operator {jfj.r^ +7i''''/i)7^ and 
hence the form factor G^{q^). This relates the previous 
results for Majorana neutrinos in Eqs. (fTT]) and 120]) . and 
clearly shows that they are indeed consistent with each 
other. 

Let us now consider matrix element calculations from 
the viewpoint of Feynman diagrams in the Dirac and Ma- 
jorana cases. An important point to note is that for 
Majorana neutrinos coupled to a graviton there are ad- 
ditional diagrams that have to be included in the cal- 
culation of {v(jp )i\J^v\v{p)i) . To be more specific we 
mean that for Majorana neutrinos there would be an ad- 
ditional "charge conjugated" diagram corresponding to 
each Feynman diagram of the Dirac case. Based on this 
observation we can try to understand how the matrix el- 
ements in the Majorana case may be calculated from the 
Dirac neutrino calculation. Again, we modify an argu- 
ment used in the neutrino- photon case Q . Let us assume 
that for Dirac neutrinos the t-channel matrix element at 
one-loop is 

{v,{p ,s)\/\J^p\ v,{p,s)),_,^^^ = K T^f3ig,j5,P,q) (29) 

where q ^ p — p and g is the electrowcak coupling 
constant. For the Majorana case the charge conju- 
gate diagram will have all particles replaced by their 
charge conjugates and the clcctroweak projection opera- 
tor Pl — *■ Pr- Unlike in the electromagnetic case (with 
coupling constant e) the coupling constant in the present 
case, namely k, does not change sign under charge con- 
jugation and hence for the corresponding Majorana case 
we have at one-loop 

(Mp ,s )\^Jafi\ '^tiP,s))uiaap = K. rf^(5,75,p,g) (30) 

+ « rf^(.9,-75,P,9) • 

This suggests that at 1-loop, compared to the Dirac neu- 
trino, the vertex factors for the Majorana neutrino will 
not have any terms proportional to 75. This conclusion 
is opposite to the equivalent conclusion in the neutrino- 
photon case [1]. 



It may be puzzling at first that this line of reasoning 
suggests in addition to the vanishing of G^^ (q^) (as in 
Eq. ([T7])) that the form factors 



Griq') - G-(g^) = . 

We would have naively expected both the argument from 
Fermi-Dirac statistics as well as the above to lead to iden- 
tical conclusions. Let us try to understand more thor- 
oughly the reason for this difference. We note that the 
1-loop argument above relied on the charge conjugation 
property of the Majorana neutrino 



= I'm 



(Majorana Condition) 



and an implicit assumption that there were no CP phases 
at 1-loop to start with in the graviton-Dirac neutrino 1- 
loop diagram. Thus it seems pertinent to explore the 
consequences of exact (or approximate) CP invariance on 
the final states, in the s-channel process, to understand 
the relation between the two arguments. 

For Dirac particles in the final state, of the s-channel 
process (iy{p )jt^(p)j|j7^iy|0), we have the well-known re- 
sult 



D /2S+1 



Lj) = (-1) 



s+i 



For clarity we have briefly reviewed the C, P and CP 
properties of Dirac and Majorana two-particle final states 
in Appendix A. Comparing the above expression with the 
angular momentum states in Eq. (jlSp we find that all 
Dirac neutrino final states are CP even except 

C ('^o) = C ('^2) = - 1 . 

For Majorana neutrino final states the CP parity is 
given by 

CA^'^'Lj) = (-1)^+1 , 

which is also the expression for the usual parity (ie. spa- 
tial reflection) since charge conjugation is trivial for the 
Majorana final states. When this is applied to the possi- 
ble final angular momentum states of the Majorana neu- 
trinos we find that all states arc CP even except 

The Majorana state is already ruled out from ar- 
guments of Fermi-Dirac statistics. So, if it is the case that 
in the Icpton sector CP invariance holds exactly, or the 
CP phases are negligible for all practical considerations, 
then from the associations in Eq. 



Gi(g2), G2{q^) : ^^o, 
G3(<z') : ^D2 



'D2 



we obtain, to all orders in perturbation theory, in addi- 
tion to Eq. dni) 



Gf( 
G\\ 



CP ^ 
CP 



G^{q') 
Gf(<Z^) 



CP 
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This is an interesting observation in the context of 
whether there is CP violation at all in the lepton sec- 
tor, especially since in the Majorana case there could 
be additional Majorana phases in the electroweak mix- 
ing matrix. One of the consequences of the above result 
is that there is a definite relation between the presence 
of CP phases (Dirac or Majorana) in the lepton sector 
and the vanishing of two of the graviton-neutrino form 
factors. 

Thus it may be speculated that irrespective of whether 
neutrinos are Dirac or Majorana, if it is feasible to probe 
the graviton form factors Gi{q^) and 6*2(9^) in future 
neutrino experiments, it could potentially yield informa- 
tion about CP phases in the charged lepton-neutrino sec- 
tor. But note that G2{q^), similar to i^2(<?^) is associated 
with a chirality flipping operator and hence suppressed 
by 0{m,y/ Ey). Thus the preliminary indications are that 
it might be challenging to probe G2{q^) in experiments, 
unless one considers very low energy neutrino sources, 
to speculate, say maybe from the cosmic neutrino back- 
ground (cf. 

Now let us consider the expressions in Eqs. ^ and 
In the Majorana off-diagonal case i ^ j, for the 
s-channel, we would have (See Appendix A) 

aC'+'Lj)=rj*Tj;i-l)^ . (31) 

Thus there are two possibilities for the s-channcl fi- 
nal state intrinsic parity, ri*ri* = — 1 and r]*ri* = +1. 
Since the product of intrinsic parities in the t-channcl 
('^(-P )j\'^tJ.'^Wip)i) is opposite to the product of intrinsic 
parities in the equivalent s-channel {iy{p )jv{p)i\J^i_i,u\0), 



the two cases of 77* ry* 



-1,-1-1 would correspond, in 



the t-channel, to the vanishing of {Ei, E2, E3) and 
(Di, D2, D3) respectively under the approximations we 
have made. This result, argued from the CP symmetries 
of the angular momentum states, is again consistent with 
the equivalent result in 

We observe specifically from Eqs. ^7}i, and ([30]) 
that the graviton vertex matrix elements for Dirac and 
Majorana cases seem to generically come out different. 
We will now show that irrespective of this apparent dif- 
ference the practical Majorana- Dirac confusion theorems 
Q first discussed in the context of neutrino scattering 
and electromagnetic interactions still hold in the grav- 
itational case and render any difference undetectable. 
The main reasons for this conclusion, as pointed out long 
back, are that the electroweak currents are left handed 
and that the neutrino masses are tiny Q . We will see in 
what follows that the arguments in the electromagnetic 
case may be essentially carried over to the gravitational 
case. 

We note that the crucial difference in the matrix ele- 
ments is that in the Majorana case we have just F^\q^) 
associated with the operator 7{p7'^} while in the Dirac 
case we have both F^{q^) and G^{q^). Thus as in the 
analogous photon vertex case if it so happens that 



d'Gfiq'- 



(32) 



to say 0{m^/E,y), then experimentally one would never 
be able to distinguish Majorana neutrinos from Dirac 
neutrinos interacting with gravity. We will in fact see 
that this is the case. Let L and R denote left and right 
chiralities. Pl and Pr will denote the electroweak left 
and right projection operators. 
At tree-level we must have 



D.O 



G: 



D.O 



(33) 



where F^'^'*^, F^'*^ and G^''^ are the terms generated at 
tree-level. This is because, crudely speaking, at tree level 
there are no additional Feynman diagrams for the Majo- 
rana neutrino compared to the Dirac neutrino. 

Now, we know that a mass term in the Dirac equation 
mixes L and R chiralities. In any electroweak loop di- 
agram there are left projection operators sitting at the 
initial and final electroweak vertices. So, in the limit of 
vanishing neutrino masses, the electroweak left projec- 
tion operators attenuate any diagram with an incoming 
or outgoing R neutrino. Thus beyond tree-level we must 
have 







(34) 



owing to the left-handed nature of the weak currents and 
the extremely small neutrino mass. The above expression 
and Eq. (|13p therefore imply 



AFfiq^) 



q'AG^iq'-) 



(35) 



where we have used the result 

Also, we have from Eq. ([30|) the observation 



AF.^' 
UL — : — UL = UL 



AFf 



q^AGi 



UL 



UL 









UL = UL 


2 


UL 



This is because, as mentioned previously, for Majorana 
neutrinos there is an additional "charge conjugated" di- 
agram, with Pl — > Pr and l~ — > Z"*", corresponding to 
each Feynman diagram of the Dirac case. 

Putting together the results (p3l) - (|36l) we have 



pD pDfl 



+9^AG^ 

pM,0 



G: 



D,0 



F" 



F. 



4 

A/,0 



T^A/.O 
_ ^3 


AF^ 


4 


+ 4 


AF^ 


AF^ 


- 4 ^ 


4 



q^AGi 



A/,0 



AFi'iq') 



Fi" 
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Thus Eq. ([32]) is indeed satisfied for the neutrino-graviton 
vertex. So, in spite of theoretical differences in the 
graviton vertex of Majorana and Dirac neutrinos, to 
0{m^l En), the two cases would be practically indistin- 
guishable in any weak-field space-time geometry. 

We note that this is in disagreement with some of the 
recent claims in the Hterature by D. Singh et. al [14| 
based on studies in the gravitational phase formalism [5| . 
Also see J. F. Nieves and P. B. Pal for a similar 

discussion that gives the answer in the negative and some 
related comments in [l5| . 

Finally, it is to be pointed out that most of our analy- 
ses for the graviton-neutrino interactions so far have been 
fairly general, based only on symmetry principles. Thus, 
we believe that most of the conclusions derived so far 
are equally applicable to other Majorana fermions that 
might exist in nature and couple to gravity. For instance, 
in the context of the proposed Minimal Supersymmetric 
Standard Model (MSSM) ^ there are additional Majo- 
rana fermions such as the gluino (g) and the neutralino 
(x°). In many models the x° plays the role of the lightest 
supersymmetric partner and may be a large component 
of the cold dark matter (CDM) in the universe. So if one 
is interested in the interaction 



X 



(37) 



in the presence of weak gravitational fields then all our 
conclusions may be carried over without almost any 
change for the Majorana neutralino 



IV. GRAVITON-NEUTRINO COUPLING AT 
1-LOOP 

We now perform a calculation of the mass off-diagonal 
neutrino-graviton form factors at 1-loop, to lowest order 
under the assumption that the neutrino energy Ei, ^ 
nil,. From this calculation we will see how the gravita- 
tional neutrino transition form factors vary with and 
also be able to estimate an i j neutrino transition 
mass radius. We are mainly interested in the mass off- 
diagonal case {i ^ j) because if it is non- vanishing it may 
be thought of as a graviton mediated "penguin" transi- 
tion Vi —> Vj for non-zero q^. Such an i — > j transition 
would be a unique quantum mechanical effect with no 
equivalent in the linearized classical theory. 

In the presence of weak gravitational fields, it is well 
known that we may use a perturbative approach to cal- 
culate the effects of the gravitational field on particles 
(see [Ol and references therein). Specifically, from the 
principle of general covariance we may derive the Feyn- 
man rules for the graviton-fermion and graviton-massive 
vector boson interactions. This is achieved by making 
the associations 



in the action integral and expanding all the relevant 
quantities around the flat metric using Eq. ^ . In 
the above, Wa are the field connections which are ex- 
pressed in terms of the vierbeins/tetrads e" 



This effective field theory approach is the one we will 
follow to study the potential effects of weak gravitational 
fields on the transition scattering of high energy neutri- 
nos. The relevant graviton Feynman rules for our calcu- 
lation are listed in Appendix B. The graviton Feynman 
rules wc derive for the standard model follow the con- 
ventions of Cheng and Li [l6j . The rules we use seem to 
be in slight disagreement with the Feynman rules used 
in (specifically regarding some of the relative signs of 
the terms), but are consistent with the rules derived by 
S.Y. Choi et. al. [l^. 



The 1-loop diagrams we are interested in for the i/i ^ Vj 
case, to lowest order, are shown in Figs. ^ and ([3]). For 
the transition form factor {Ek{q^), Dk{q'^)) calculation 
these are the only Feynman diagrams we consider. If 
we are instead interested in the diagonal form factors 
{Fk{q^), Gk{q^)) there would also be additional Feyn- 
man diagrams involving the Z boson. We choose to work 
in the Feynman- 't Hooft gauge using dimensional regu- 
larization and will adopt the M S scheme. Since we are 
working to lowest order we neglect terms that are propor- 
tional to neutrino masses, since by assumption E^, ^ m^. 
Also, we do not include diagrams with ghost fields run- 
ning in the loop which are expected to give corrections 
at a higher order. These are the approximations in our 
loop calculations. 



Our task now is to explicitly evaluate the diagrams in 
the presence of a weak gravitational field to lowest or- 
der, obtain the q^ dependence of the gravitational form 
factors and then estimate the neutrino transition mass 
radius. Let us take up the first diagram shown in Fig. 
(121). It has a topology that is unique to the graviton 
interaction and is not present in the photon case. As 
before we denote the neutrino mass-eigenstates by the 
latin indices i,j and I denotes the lepton flavor. The 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [ll| 
is denoted by Uu. The mass and the mass of the 
charged lepton l^ are denoted by Mw and m; respec- 
tively, g is the electroweak coupling constant. 



Using the Feynman rules for linearized gravity, as re- 
viewed in Appendix B, the spin-2 gauge current for the 
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first Fcynman diagram with i ^ j in Fig. ([2|) is 

df^k _ , „ , ig 
■V2 



E 



, ^ + mi ^ , Kg X \ / \ 



2^/2 



,M2 



ng^U^jUu 



Uj{p') 



(27r)4 (fc2 -m2)((fc-p')^ 



Ml) 



where 



For the similar Feynman diagram in Fig. ([2|) , where the 
graviton is coupled to the final vertex we have 



E 

l—e.fi,T 



Uj{p') 



4 

(27r)4 {k^ ^ mf){{k - py 



Ml) 



The above spin-2 gauge currents for the diagrams in 
Fig. (12) may be simplified using the standard techniques 
of Dirac algebra, assuming E^, 3> mjy. This gives finally 
the approximate expression for the gauge current 



^ 327r2 



Ul.Uu / dx (38) 



where 



A{Ml,mf) = xMl + (1 - x)mf 



Studying the above 1-loop expression we see that in- 
deed we do get an operator form that is consistent with 
our expectations of gauge invariance. The only form 
factors that are generated at this order are E-i{q^) and 
D3{q^). It is noted that none of the other form factors 
are generated at this order from these diagrams. As pre- 
viously mentioned, the form factors of the type E2{q^) 
and D2{q^) corresponding to TaVp mediate chirality flips 
for the incoming neutrino and are thereby further sup- 
pressed by 0{m^ / Ej^). Since our approximate calcula- 
tion neglects terms proportional to nii, these terms are 
not retained. 




d k _ f.. Kg Xr> TT* 



FIG. 2: Four-point neutrino graviton diagrams at 1-loop with 



Also note that the contributions from the four-point 
graviton diagrams do not have a q^ dependence at this 
order of approximation. As we shall see, this implies that 
when we enforce correct renormalization for the spin-2 
currents they would not contribute. Thus this would im- 
ply that they do not contribute to the vi Vj transition 
mass radius of the neutrino as well at leading order. 

The last thing we would like to point out is that 
the Glashow-Iliopoulos-Maiani (GIM) mechanism [l^l is 
again relevant here. Note from the above that in the limit 
of degenerate charged lepton masses, the loop result van- 
ishes for i ^ j ■ 

The next diagram we consider is the one where the 
graviton is coupled to the gauge boson in the loop. 
In the equivalent photon case, such a diagram would only 
have contributed to next order in 1/M^. A quick look at 
the Feynman rule for the gauge boson-graviton coupling 
would convince us that this is not true in the graviton 
case. This is because in the graviton case the numerator 
of the gauge boson-graviton coupling has a term that is 
proportional to Myy. 

For the graviton coupling to the W boson we have 



E 



d^k ^ig^u;,UH 

(27r)4 2 



(p + fc)2 _ Mlip + k'f - Mlik^ - mf) 



11 



where 



+ p'pP{y-nt,}a-'p'{i,Py}-nai3~ {p' ■p)-nuia'ny}i3\ ■ 

Once again, the expression above may be simphfied 
using the standard techniques of 1-loop calculations. For 
brevity, again we present only the final result. The W- 
boson graviton coupling of Fig. ^ gives a contribution 
to the gauge current 



AJi™*)u,^E^f/yC/H / dxdy (39) 



(3-2x-2j/)log [A'(g2,M, 



xyq 



T-7±(p'+k) 





r{p-+k) 
r{p+k) 



where 



FIG. 3: Graviton-W boson and graviton-charged lepton cou- 
pled diagrams at 1-loop for the case i ^ j- 



A {q\ M^,mf ) = (x + y)M^ + {1 - x - y)m: 



xyq 



In the integral, R denotes the region of the X — Y plane 
bounded by the axes and the line x + y = 1. Again, we 
note that the rank-2 tensor current comes out to have the 
required gauge invariant form. Note that the contribu- 
tion from this diagram does have explicit q^ dependences 
at this order. Hence we would expect it to contribute to 
the q^ dependence of E'i[q^) and the neutrino transition 
mass radius. The GIM mechanism [l^ again suppresses 
the contribution for i ^ j. 

The last Fcynman diagram wc consider is the one 
where the graviton is coupled to the charged lepton in 
the loop, as shown in Fig. A similar diagram also 
gives a contribution in the photon-neutrino coupling case. 
For the diagram under consideration the tensor current 
is then given by 



l—e,fj,,T 



2 ./ {2ttY ' 



2\ « 



(P) 



{{p' + ky - m^){k^ - M^){{p + fc)2 - mf) 
where 

= ^[i{t.{p + p'),^}]-\vt^Ap + f-'^^'^f] ■ 

The above expression may be grouped in powers of mi 
and simplified. It is found that the only terms that sur- 
vive arc those proportional to mf and constant terms. 



Simplifying the expressions and performing the momen- 
tum integration finally gives 



[x + y-l) 



-ilog [^\q\M^,mt)\ 
{l-x~y + xy)q 



(40) 



21 



2A"{q^,M^,mf) A"iq^,M^,mf) 
where 

A"{q^,M^,mf) = {x + y)mf + {l-x-y)M^-xyq^ 

Again in the Fcynman integral, R denotes the region of 
the X — y plane bounded by the axes and the line x + y ^ 
1. 

Now we may put together the contributions from all 
the Feynman diagrams to get the net spin-2 tensor cur- 
rent. The combined 1-loop gauge current is 

AJ^V°°P(g^) - AJ^r-iq^) + AJ^,* (g2) + A (g2) 

(41) 

where AJ[q^) denotes the tensor current renormalized 
according to 

AJ^,(g2) = AV(g2) _ AV(0) 



in accordance with the conditions in Eq. pop . Thus we 
observe that the contributions from the four-point gravi- 
ton coupling in Fig. ^ is renormalized away since it has 
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FIG. 4: \E3{q^)\ Vs. -q^ plot for i / j illustrating the ap- 
proximate dependence to of the dominant transition form 
factor at lowest order. Note that the q^ dependence is highly 
suppressed and we observe significant deviations from the van- 
ishing q —* value only at very large values of momentum 
transfer. 



no dependence at the order we are considering. The 
dominant contribution to E^{q^), D^{q^) and hence to 
the i ^ j transition mass radius, at this order of approx- 
imation, therefore comes from the Feynman diagrams in 
Fig. ([31). Also, as mentioned before, the off-diagonal 
contribution is suppressed as should be expected by the 
GIM mechanism [l^ . In the limit of degenerate charged 
lepton masses the 1-loop contributions would completely 
vanish and there would be no i ^ j transitions. 

In Fig.Q we plot how, to leading order, the form factor 
^3(9^) [* ^ j] changes with momentum transfer ~q^. It 
is seen that the q^ dependence for moderate values of 
momentum transfer is very weak and almost insignificant. 
If there is any possibility of checking this effect in the 
gravitational scattering of neutrinos, one would have to 
consider some effect where the momentum transfer is very 
far from 0. For the scenarios we consider where ^ 
TO,,, note that 



q^ Ki - 4El sin^ 



Thus larger and larger values of momentum transfer cor- 
respond to larger and larger scattering angles 9 for fixed 
values of Ejj. 

Let us now estimate a neutrino transition mass radius 
based on Eq. (|1T|) . Crudely speaking, the diagonal mass 
radius (in mass basis) of a neutrino may be defined as the 
second moment of a hypothetical neutrino mass density 



cPr \Tf p"^{r) 



(42) 



Here as in section I, p"^{r) is the fictitious neutrino mass 
density. If we assume that the hypothetical neutrino 
mass-density is spherically symmetric and that the mo- 
mentum exchange q is small then we may re-write the 



above expression as in Eq. (fTT|l . This result may now be 
assumed to hold true even when Eq. (|42|) is not strictly 
correct and when i ^ j- This is again in parallel with the 
arguments leading to the definition of a neutrino charge 
radius. In our case then, a — > Vj transition mass radius 
may be defined as 



E2{q^) is not generated in our leading order calcula- 
tion and E'i{q^) may be calculated from the expressions 
we obtain, after numerically evaluating the Feynman in- 
tegrals. Using Eqs. ((38l) - (|4T|l to calculate the above, 
say for V2 — > I'a, we find the approximate value for the 
Vi Vj transition mass radius 

(r2)^^_^^j2 ^ 3] - 1.5 X 10-3^ cm^ 

A neutrino transition charge radius calculation with 
the same assumptions as the one we adopt (again for 
2^3) gives 



1.1 X 10"^^ cm^ 



Now, it is known that the neutrino charge radius is 
gauge dependent and hence unphysical, but there have 
been attempts to arrive at a gauge invariant definition 
for the neutrino charge radius. (See for example the rele- 
vant section in [l^l and references therein.) Also, all the 
best-fit values quoted for the neutrino charge radius are 
in the ball-park of 10^^^ — 10^^^ cm^ [l^l- Our approxi- 
mate calculation of a neutrino transition mass radius thus 
seems to indicate that 



< 



There is an intuitive way to understand the difference 
in order of magnitudes between the neutrino transition 
charge radius and neutrino transition mass radius. Let 
us consider the vacuum polarization of the neutrino. For 
the neutrino transition charge radius, the photon couples 
to the vacuum polarized and with equal "weight" , 
since they have the same magnitude of charge ±1. Now, 
the definition of the Compton wavelength is 

A. A 

TOC 

where h is the Planck's constant, c is the speed of light 
and TO is the mass of the particle. Thus the charged 
lepton has a larger Compton wavelength compared to the 
boson and thereby it contributes more significantly 
to the neutrino charge radius. Taking a charge weighted 
average, we may associate heuristically 



19/- 



|Af- + k^+|A2 



(43) 



where qi- and qw^ denote the respective charges and Aj- 
and \w+ a-rc the Compton wavelengths of l~ and . 
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Similarly, in the transition mass radius probed by the 
graviton, taking a weighted average with the respective 
masses, we may claim that 



w+ 



mi 



(44) 



approximately. 

Calculating the Compton wavelength for the charged 
leptons (appropriately weighted by the relevant PMNS 
matrix factors for the 2^3 case) and the W±, we get 
from Eqs. (gSl) and (|44]) 



10" 



Thus, it is seen that the value calculated from heuris- 
tic arguments is in the ball park of the equivalent ratio 
calculated from 1-loop calculations, 



(ri): 



\ I' /charge 



10" 



This probably gives a conceptual way of understanding 
the smallness of the neutrino mass radius compared to 
the neutrino charge radius. 



V. SUMMARY 

In this paper we studied the interaction of neutri- 
nos with gravitational fields in the weak field regime. 
Among our aims were to understand the symmetries of 
the graviton-neutrino vertex in some detail and explore 
the gravitational neutrino transition form factors gener- 
ated at 1-loop. 

The form of the graviton-neutrino vertex was moti- 
vated based on Lorentz and gauge invariance for both 
the mass diagonal and off-diagonal cases. The non- 
relativistic interpretations for the neutrino gravitational 
form factors were derived. Also, renormalization crite- 
ria on the matrix elements in the mass basis were im- 
posed so as to preserve the weak equivalence principle. 
Wc then used general arguments of spin-statistics, CP in- 
variance and most importantly the symmetries of the an- 
gular momentum states to deduce theoretical differences 
between the Majorana and Dirac cases. Specific associa- 
tions were made between the neutrino gravitational form 
factors and angular momentum states that arc allowed. 
Differences in angular momentum states for the graviton 
case compared to the photon case were pointed out, for 
instance, the presence of states due to the J — part of 
the off-shell graviton. 

It was then proved that, as in the electromagnetic case, 
the practical Majorana-Dirac confusion theorems are still 
valid in the gravitational case. This meant that in spite 
of the theoretical differences, due to the tiny neutrino 



masses and left-handed nature of the electroweak cur- 
rents, the two cases (Majorana and Dirac) would be ef- 
fectively indistinguishable for any space-time geometry 
satisfying the weak field condition. Finally we made some 
brief remarks about the results in sections II and III being 
valid for other Majorana particles that may be present 
in nature. 

We then calculated the neutrino-graviton transition 
form factors generated at 1-loop to leading order and ob- 
tained an approximate understanding of their depen- 
dence. We pointed out that the gravity mediated i ^ j 
effect would be purely quantum mechanical and absent in 
the classical theory of linearized gravity. A neutrino tran- 
sition mass radius in the presence of a gravitational field 
was then estimated from the approximate 1-loop expres- 
sions and compared to the neutrino transition charge ra- 
dius. It was found from our calculation that the neutrino 
transition mass radius generally comes out to be smaller 
than the neutrino transition charge radius by couple of 
orders of magnitude. Finally, we tried to give a physical 
explanation for the smallness of the neutrino transition 
mass radius compared to the neutrino transition charge 
radius. 

There are many avenues left to be explored relating to 
the study of high energy neutrinos interacting with grav- 
ity. A more detailed calculation of the neutrino gravi- 
tational form factors with higher order terms included 
should clarify whether the behavior seen in our approxi- 
mate computations for the neutrino transition form fac- 
tors and neutrino transition mass radius is robust. Re- 
lated to this is the question of gauge invariance of the 
gravitational form factors and whether they may even in 
principle be directly measured in some future neutrino 
experiment. This may be relevant in the context of ex- 
ploring CP phases in the lepton sector as we remarked 
in section III. Other speculations are whether we may 
use the differences in gravitational form factors of parti- 
cles as cosmological probes, say for example to ascertain 
whether dark matter is bosonic or fermionic. 
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APPENDIX A: PARITY AND CHARGE 
CONJUGATION PROPERTIES FOR DIRAC AND 
MAJORANA TWO-PARTICLE FINAL STATES. 

For completeness we briefly derive the well known P, 
C and CP properties of two particle fermion states in the 
Dirac and Majorana cases. 



2. Majorana two-particle final states 

The bound state of the Majorana particles may be rep- 
resented to lowest order again as 

\vv-P,S)m = 51 Z*^'^ <^M{k-si,S2) 4,,,a^^j0) 



1. Dirac two-particle final states 

For the bound state of two Dirae particles we may 
represent the final state to lowest order as 

\vD;P,S)D=y^ I <i'k^D{k;si,S2) ai |0) + . . . 

Sl,S2 

Let us first consider the effect of the parity operator V 
on the final state. Using the fact that (see for example 

M) 

Vai V = 111 r 

VP - V = ^I- 

and 77ft = — ry* we have after re-defining the variable of 
integration 



In the above expression 

$£,(fc ;si,S2) = -^D{-k;si,S2) 

Using the familiar properties of the spherical harmonics 
this gives the well known result 

Here L and S denote the angular momentum quantum 
number and total spin. J is the total angular momentum 
defined in the usual way as J = L + S. 

Similarly, applying the charge conjugation operator C 
on the Dirac final state gives 

C\i^D-P,S)d ^Y. J " '^'oik" -,81,82) 4",.i^L-"..jO> 



where 

^jj{k ;si,S2) = -*n(-^;s2,si) 

Using the properties of spherical harmonics and the spin 
states we get the familiar result 

From the above two results it is easy to see that the 
CP parity for the Dirac two-particle final state is given 

by 

C C'^'Lj) = (-1)^+^ 



where we have implicitly used the Majorana condition. 
Acting on the state by the parity operator V now gives 



Unlike Dirac particles, Majorana particles have imagi- 
nary intrinsic parity. This is most easily seen by consid- 
ering the expression M{x,t)V~^ and requiring it to 
be proportional to ^'7v/(— a?, t). This gives the condition 



-ila 



implying that rja = ±i. Using this condition, after a 
re-definition of the momentum integration, we get 

V\vv-P,S)m [ (i'k $;,(fc';si,S2) 4' «^ r |0) 



where 

<^M{k ; si, S2) = si, S2) 

which is similar to the Dirac case. This gives again 

The action of the charge conjugation operator on the 
Majorana final state does not yield any new information 
since we have 



Thus 



C\vv-P,S)m = \vv;P,S) 



Q C^+'Lj) = +1 



M 



Putting the above results together we get finally for 
the CP parity of Majorana two-particle final states 



c 



M (2S+1 



Lj) = (-1) 



L + l 



If the Majorana particles in the final state were two 
different eigenstates a ^ h then in the above expressions 
VlVb 7^ ^1 ^^"^ parity would explicitly depend 

on the intrinsic parity of each particle. Thus in the off- 
diagonal case we would have 

ar+iL,) = ^:%*(-i)^ 

From the above results we see that there are some key 
differences in the Majorana two-particle final state com- 
pared to the Dirac two-particle final state for charge con- 
jugation, parity and CP. 
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APPENDIX B: RELEVANT FEYNMAN RULES 
DERIVED FOR LINEARIZED GRAVITY 





ft'-(q) 



FIG. 7: Graviton-fermion coupling is — |^[7{fi( 



FIG. 5: For the four point graviton coupling the Feynman 
rule is i0= [ri^„riai3 - ^ri^fjri„a - ^Tl^aTl^isjY^ Pl 




FIG. 6: The Feynman rule is it— ^ [?7fii^7?a,3 — Vt^arivp - 

?7^a')M/3] -*f hfif (P' -P VafI -Po^P'p) +P'{^,PaV0-^} +P'l3P{<^VlJ.}a ^ 

p'(^,Py}riafi - {P •p)'7{mq»7>'}/3] 
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